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Abstract. 

The fourth harmonic of the azimuthal distribution of particles t>4 has been measured 
for Au-Au collisions at the Relativistic Heavy Ion Collider (RHIC). The centrality 
dependence of V4 does not agree with the prediction from hydrodynamics. In particular, 
the ratio V4/(v2) 2 , where «2 denotes the second harmonic of the azimuthal distribution 
of particles, is significantly larger than predicted by hydrodynamics. We argue that 
this discrepancy is mostly due to elliptic flow (i^) fluctuations. We evaluate these 
fluctuations on the basis of a Monte Carlo Glauber calculation. The effect of deviations 
from local thermal equilibrium is also studied, but appears to be only a small correction. 
Combining these two effects allows us to reproduce experimental data for peripheral 
and midcentral collisions. However, we are unable to explain the large magnitude of 
V4/(v2) 2 observed for the most central collisions. 
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1. Introduction 



The azimuthal distribution of emitted particles is a good tool for understanding the 
bulk properties of the matter created in non central nucleus-nucleus collisions. In the 
center of mass rapidity region, it can be expanded in Fourier series: 
dN 

— oc 1 + 2v 2 cos(20) + 2v A cos(40) H (1) 

where is the azimuthal angle with respect to the direction of the impact parameter, and 
odd harmonics are zero by symmetry. The large magnitude of elliptic flow v-i observed 
at RHIC suggests that the matter created in Au-Au collisions behaves like an almost 
perfect fluid. However, recent experiments [H[2] observe that, at midrapidity and fixed 
p t , t> 4 ~ (t>2) 2 , while ideal hydrodynamics predicts that V4 = \{v2) 2 [3]- I n this talk, I 
investigate this discrepancy. 

2. Fluctuations in initial conditions 




Figure 1. (Color online) Picture of the two frames used for denning the initial 
eccentricity (from [4]). The x axis defines the reaction plane while the x' axis is 
the minor axis of the ellipse drawn by the participating nucleons (grey dots). 

Figured] presents a schematic picture of a non central heavy-ion collision (HIC). The 
overlap area of the nuclei has an almond shape, which generates elliptic flow. However, 
the matter is not continuously distributed in a nucleus. The positions of the nucleons 
in the colliding nucleus are important: they also draw an ellipse which differs from 
the overlap area both in eccentricity and in orientation. From one event to the other, 
even at fixed impact parameter, the positions of the nucleons in the nucleus fluctuate. 
The participant plane eccentricity (epp), defined as the eccentricity of the ellipse drawn 
by the participating nucleons [SI E], thus fluctuates. Since elliptic flow appears to be 
driven by this participant plane eccentricity, these eccentricity fluctuations translate 
into fluctuations of the flow coefficients v 2 and V4 [7]. 
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3. Modeling eccentricity fluctuations 



The initial distribution of energy, which is needed to compute the initial eccentricity in 
a HIC, is poorly known. In this talk I use a specific model, based on a Monte Carlo 
Glauber (MCG) calculation [8]. The initial eccentricity is given for each event by: 



v K 2 -^) 2 + 4< 

€PP = ^ 

where a 2 = (x 2 ) — (x) 2 and a xy = (xy) — (x)(y) and the () denote averages over 
participating nucleons. Each participant nucleon is given a weight proportional to 
the number of particles it creates, according to the two-component picture: w = 
(1 - x) + xN coU ^ nudeon where N coU _ nudeon is the number of binary collisions of the 
nucleon. The sum of weights scales like the multiplicity: 
dN ch 

= n„„ 11 — t. 1 — 

2 



x) —P^l + xN co u 



drj '" w r ~ J 2 ' ~" CU "J ' (3) 
where N part and N coU are respectively the number of participants and of binary collisions 
of the considered event. We choose the value x = 0.13 which best describes the charged 
hadron multiplicity observed experimentally [9]. We define the centrality according to 
the number of participants. We evaluate eccentricity fluctuations in centrality classes 
containing 5% of the total number of events. We do not introduce any hard core 
repulsion between nucleons in the MCG. 



4. How eccentricity fluctuations affect t> 4 /t> 
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Figure 2. (Color online) Centrality dependence of Vi/(v2) 2 '- data from STAR [10] and 
PHENIX jTT] ; error bars on STAR data points are our estimates of nonflow errors [12] • 
Lines are predictions from ideal hydro with or without fluctuations. 

There is no direct measure of the flow coefficients V2 and They can be 
obtained using different analysis methods. The one I will consider now relies on 
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azimuthal correlations between particles near midrapidity. Experimentally, t> 2 can be 
extracted from the 2-particle correlation and V4 from the 3-particle correlation using 
(cos(20i — 20 2 )) = ((^2) 2 ) and (cos(40i — 20 2 — 203)) = (f4(f2) 2 ), where angular brackets 
denote an average value within a centrality class. Thus, any experimental measure of 
V4/V2 obtained using this method is rather a measure of (1*4)/ (v 2 ) 2 . Taking into account 
the ideal hydrodynamics prediction t> 4 = |(v 2 ) 2 [3], we obtain: 

( v 4 \ 1 ((v 2 )*) > 1 



measured 

2 (fa)*)* 2 

Assuming that v% scales like the participant plane eccentricity epp, the effects of 
fluctuations on v^jv\ is obtained by computing: 

( \ 1 (4p) 

\WJmoo ~ 2 ^pp) 2 ' 
The resulting prediction for v^jv\ is displayed in figure [2J Fluctuations clearly 
explain most of the difference between hydro and data. It also appears that experimental 
data are still slightly higher than our prediction from fluctuations. However, these 
predictions are based on a specific parametrization of the initial conditions. 

5. Flow fluctuations from experimental data 
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Figure 3. Same as figure [21 additional points (labeled "fluct. from V2") are obtained 
using equation ([6|) with W2{4} from [13] and i>2{2} from [14]. 

Another possible way of evaluating flow fluctuations is to compare the values of v 2 
obtained using different analysis methods. Elliptic flow can be obtained from both 2- 
particle cumulants (w 2 {2}) and from 4-particle cumulants (t> 2 {4}) using w 2 {2} 2 = ((w 2 ) 2 } 
(neglecting the non-flow contribution) and f 2 {4} 4 = 2((f 2 ) 2 ) 2 — ((f 2 ) 4 ). Inverting the 
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last equation leads to: 

\ _ 1 ((V2) 4 ) 1 / /^{4}\ 4 \ 

"2((^ 2 ) 2 ) 2 2 ^ \v 2 {2}) J' 1 ) 

The values of t^/vf obtained using this method are displayed on figure [3J They 
overshoot slightly our results from MCG eccentricity fluctuations, but the overall 
agreement remains good. This provides a good check of our MCG prediction. A small 
residual discrepancy remains between our prediction and the experimental data. We 
argue that for peripheral to midcentral collisions, it may be understood in terms of 
deviations from local thermal equilibrium. 



6. Partial thermalization effects 



So far, I have only discussed how fluctuations in initial conditions modify the prediction 
from ideal hydrodynamics for v±jv\. But ideal hydrodynamics relies on the very strong 
assumption that the system remains in local thermal equilibrium (a regime where the 
average number of collisions per particle n co u is large) throughout the evolution. In a 
previous work [T5] we have shown that, in order to reproduce the centrality dependence 
of elliptic flow, the deviation from local thermal equilibrium must be taken into account 
{n co u oc 3 — 5 would be a typical value for Au-Au collisions at the top RHIC energy). 



11 

.. • ■ - 

transport results 1 — ■ — 1 
linear fit 

0.2 0.4 0.6 0.8 1 1.2 
K 

Figure 4. (Color online) Variation of V4/(v2) 2 with the Knudsen number. 
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Qualitatively, in the limit of small n co u (far from equilibrium), one expects both 
t> 2 and t> 4 to scale like n co u, so that v^jiyq) 2 scales like l/n coU : we thus expect that the 
farther the system from equilibrium, the larger Vi/^2) 2 [16j. In order to have a more 
quantitative estimate of the effects of partial thermalization, we use a 2 + 1-dimensional 
solution of the relativistic Boltzmann equation to study systems with arbitrary n co u. We 
use the Knudsen number [16], K oc l/n co u, as a measure of the degree of thermalization 
of the system. 
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Figure H] displays the dependence of v^/^) 2 with the Knudsen number. In the 
limit K — > 0, transport results show that v^/^) 2 = 0.52, which is close to 1/2. We also 
observe, as expected from the low n co u limit, that increasing K leads to an increase of 
Vi/{v2) 2 - But this effect is only a small correction. 
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Figure 5. Same plot as figure [2] (right), with one additional curve showing the effect 
of the deviation from local equilibrium. 

The effects of partial thermalization on the centrality dependence of v<i/(v2) 2 are 
displayed on figure [51 The values of the Knudsen number needed for this plot are 
borrowed from a previous study [15]. Figure [5] shows that adding the effects of deviation 
from local thermal equilibrium to the fluctuations, our prediction overshoots slightly the 
data for midcentral and peripheral collisions, but the overall agreement is good. We do 
not yet understand the large value of Vi/{v2) 2 for central collisions. 

7. Conclusion 

To conclude, I would like to recall three points: 1) V4 is mainly induced by V2] 2) 
the deviation from local equilibrium has a small effect on f4/(f2) 2 ; 3) eccentricity 
fluctuations explain the observed values of v^j (v 2 ) 2 , except for the most central collisions 
which require further investigation. 
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